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ABSTRACT 

Aims. In this work we study time-dependent precipitation of an electron beam injected into a flaring atmosphere with a converging 
magnetic field by considering collisional and Ohmic losses with anisotropic scattering and pitch angle diffusion. Two injection regimes 
are investigated: short impulse and stationary injection. The effects of converging magnetic fields with different spatial profiles are 
compared and the energy deposition produced by the precipitating electrons at different depths and regimes is calculated. 
Methods. The time dependent Fokker-Planck equation for electron distribution in depth, energy and pitch angle was solved numeri- 
cally by using the summary approximation method. 

Results. It was found that steady state injection is established for beam electrons at 0.07-0.2 seconds after the injection onset depend- 
ing on the initial beam parameters. Energy deposition by a stationary beam is strongly dependent on a self-induced electric field but 
less on a magnetic field convergence. Energy depositions by short electron impulses are found to be insensitive to the self-induced 
electric field but are strongly affected by a magnetic convergence. Short beam impulses are shown to produce sharp asymmetric hard 
X-ray bursts within a millisecond timescale often observed in solar flares. 
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1. Introduction 

Observations of solar flares in hard X-rays provide vital in- 
formation about scenarios, in which accelerated electrons gain 
and deposit their energy into flaring atmospheres. In recent 
years the theory describing the generation of bremsstrahlung 
emission has been significantly progressed in many direc- 
tions by improving the mechanisms for emitting this radia- 
tion, e.g. considerin g relativistic bremsstrahlung cross-sections 
dKontar et alJ I2006I) . taking into account various aspects of 
the photospheric albedo effects while deriving mean elec- 
tron spe£tra_Jrcm die observed bremsstrahlung photon spec- 
tra dKontar et alJ I2006I) . On the other hand, substantial im- 
provements were also achieved in the solutions of the direct 
problem of electron precipitation into a flaring atmosphere by 
taking into account differen t mechanisms of electron energy 
losses: Coulomb collisions dBrownl Il97lb iBrown et all [2000) 
combined with the deceleration by the self-induced ele ctric field 
( Zhar kova etaf] 1 1995b Izharkova & Gordovskyyl [2006) in flar- 
ing atmospheres with strong temperature a nd density gradients 
derived from the hydro-dynamic solution s ( Somov et al.l 1981, 
ll982HNagai & Emslid[T984L iFisher et al.lll985l) . 

Substantial progress in the quantitative interpretation of hard 
X-ray emission has been made in recent years by using high 
temporal and spatia l resolution ob servations carried out by the 
RHESSI payload dLin et al.ll20"03h . The latter provides the lo- 
cations and shapes of hard X-ray sources on the solar disk, 
their temporal variations and energy spectra evolution du ring the 
flare duration dHolman et alJl2003t iKrucker et al 1 l2008h . These 
observations are often accompanied by other observations (in 
microwaves (MW), EUV and optical ranges) which reveal a 
very close temporal correlation between H XR and MW, UV 
and even optical em i ssion (see for exampl e, iKundu et al.ll2004t 
iFletcher et al"1l2007t iGrechnev et a l. 2008). This highlighted a 



further need for the improvements of electron transport mod- 
els, which can simultaneously account both temporarily and spa- 
tially for all these types of emission. 

The RHESSI observations of double power la w energy spec- 
tra wi th flattening towards lower photon energies (Holma net al.l 
120031) . which leads to the soft-hard-soft temporal pattern o f 
the photon spectra indices below 35 keV (Grigis & Benz 2006), 
highl ighted their role in the formation of the self-induced electric 
field (Zh arkova & Gordo vskvv 2006). The authors considered a 
stationary beam injection and naturally reproduced such a spec- 
tral flattening by electron deceleration in the electric field, in- 
duced by beam electrons themselves. The flattening was shown 
to be proportional to the initial energy flux of beam ele ctrons 
and their spectral indices (IZharkova & Gordovskyyl20 06). Then 
the soft-hard-soft pattern in photon spectra above can be easily 
reproduced by a triangular increase and decrease of the beam en- 
ergy flux in the time interv al of a few seco nds, which was often 
observed both by RHESSI ( Lin et al. 2003), and by the previous 
SMM mission dKane et al.ll980h . Furthermore, numerous obser- 
vations of solar flares by SMM, TRACE and RHESSI suggest 
that the areas of flaring loops decrease and, thus, thei r magnetic 
fields increase with the depth of the solar atmosph ere dLang et al.l 
119931: iBrosius & Whitdl2006l: iKontar et al.ll2008l) . This increase 
of magnetic field can act as a magnetic mirror for the precipi- 
tating electrons forcing them to return back to the source in the 
corona, in addition to self-induced electric field. 

In the present paper we propose two models of the magnetic 
field variations with depth. One model is a fi tting to the mea- 
surem ents of the magneti c field in the coron a dBrosius & Whitd 
2006) and chromosphere dKontar et al.ll2008l) . another shows the 
exponential increase of magnetic field from the corona to the 
upper chromosphere while remaining a constant in the lower 
chromosphere. The outcome is compared with the two other 



2 



T. V. Siversky & V. V. Zharkova: Stationary and impulsive injection of electron beams 



models proposed earlier: the first one dLeach & P etrosian 1981) 
where the authors assumed that the magnetic column depth scale 
<9 In B/ds is constant implying the exponential mag netic field in- 
crease with a column depth and the second one (iMcClementsI 
119921) considering a parabolic increase of the magnetic field with 
a linear depth. 

Also the temporal intervals of impulsive increases of 
HXR emission vary from very short (te ns of milliseconds 
dKiplinger et al.|[l98 3; Chari kov et al.ll2004l)) to t ens of minutes 
often observed by RHESSI (iHolman et al.ll2003h . This encour- 



aged us to revise the electron transport models and to consider 
solutions of a time-dependent Fokker-Planck equation for dif- 
ferent timescales of beam injection (milliseconds, seconds and 
minutes). The electron transport, in turn, can slow down also 
by anisotropic scattering of beam electrons in this self-induced 
electric field enhanced by their magnetic mirroring in converging 
magnetic loops. The further delay can be caused by the particle 
diffusion in pitch angles and energy which can significantly ex- 
tend the electron transport time into deeper atmospheric layers 
where they are fully thermalised. 

We also apply the time dependent Fokker-Planck equation in 
order to compare the solutions for electron precipitation for sta- 
tionary and impulsive injection and their effect on resulting hard 
X-ray emission and ambient plasma heating for different param- 
eters of beam electrons. We also investigate these Fokker-Planck 
solutions for the different models of a converging magnetic field 
by taking into account all the mechanisms of energy loss (colli- 
sions, Ohmic losses) and anisotropic scattering but without dif- 
fusion in energy. 

The problem is formulated in Sect.|2]and the method of solu- 
tion is described in Sect. [3] The stationary injection into flaring 
atmosphere with different magnetic field convergence and colli- 
sional plus Ohmic losses with anisotropic scattering is consid- 
ered in Sect. [4] and the impulsive injection for short timescales 
below tens milliseconds is considered in Sect. [5] The discussion 
and conclusions are drawn in Sect. [6] 



2. Problem formulation 

2.1. The Fokker-Planck equation 

We consider a one-dimensional beam of high energy electrons, 
that is injected into the solar atmosphere. The beam electron 
velocity distribution /, as a function of time t, depth /, veloc- 
ity V and pitch angle between the velocity and the magnetic 
field 6, can be found by solving the Fokker-P lanck equation 
dDiakonov & Somovl 198 ft IZharkovaet ai1ll995l) : 
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where the collisional rate v is given by 
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£ is the self-induced electric field, B is the background mag- 
netic field, n p (l) is the density of the ambient plasma, In A is the 
Coulomb logarithm, e and m e are the electron charge and mass 



respectively. In our study we assume that the Coulomb logarithm 
is constant: In A m 20. 

Let us introduce the following dimensionless variables: 
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where £ is the column depth: 
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Eq = 12 keV is the lower cut-off energy and no = 10 cm 
Eq. dTJ in dimensionless variables takes the form 
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where a% is the magnetic convergence parameter defined as 
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As it was shown by Ivan den Oordl d 19901) after a few colli- 
sional times from the start of the beam injection the background 
plasma fully compensates the charge and current deposited by 
the beam. In our model we assume that the beam current is al- 
ways compensated, i.e. jt, + j p = 0, where y'j and j p are the 
beam and plasma currents respectively. Note, that beam current 
jb includes electrons with ft < 0, i.e. those that moves back to 
source, which can significantly decrease the compensation cur- 
rent j p . Then, by applying the Ohm's law for j p and taking into 
account that j„ = — jt, we can find the self-induced electric field 
(iMcClement s 1992). 



e (t, s) 



(T 



^min -1 



where the dimensionless conductivity is 

o- (s) = 1 .97 V27T- — — 

4 ymeFo 



3/2 



(12) 



(13) 



and Fo = 10 10 erg crrT 2 s~' is the normalisation factor for the 
energy flux of the beam. The ambient plasma is assumed to be 
preheated and its density, n, and temperature, T, as functions 
of the column depth, s, are calculated using the hydro-dynamic 
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is equal to some preset value F top , where Fq - 10 lu erg cm 2 s 1 
is the normalisation factor of the energy flux. At large depth, 
s = s max = 9.17 x 10 2 (or 1.01 x 10 23 cnr 2 ), the number of 
electrons in the beam is assumed to be negligibly small, thus the 
corresponding boundary condition is 
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Fig. 1: Density and temperature of the ambient plasma calcu- 
lated by the hydro-dynamic model (Zharkova & Zharkov 2007) 
for different beam energy flux, Fq, and power law index, y, of 
the beam electron distribution (see Eq. $15[). 



model (Zharkova & Zharkovll2007h . The profiles n(s) and T(s) 
for different beam parameters are shown in Fig.Q] These profiles 
do not change in time, since the thermal conduction processes 
have much longer time scales than the precipitation processes 
studied here. 

The electron beam as well as the plasma return current 
can be a subject of vari ous plasma instabilities. For example, 
Zhar kova & Ko bvlinskii (1 1992b have shown that a beam with 
energy flux higher than 10 11 erg cirT 2 s" 1 can be unstable with 
respect to the ion-sound waves. However, such instabilities are 
out of scope of the current paper and will be considered in forth- 
coming studies. 

2.2. Initial and boundary conditions 

There are no beam electrons before the injection starts, thus, the 
initial condition is: 



/(r = 0, s,z,fi) = 0. 



(14) 



The boundary condition at s = s m j n = 2.08 x 10 3 (or 
2.29 x 10 17 cm" 2 ) corresponds to the injected beam distribution 



,(5-1 



f(T,S= S min ,Z,fi > 0) = fnlff(T) 



z s+y + l 



exp - 



A// 2 



(15) 

where A/u is the initial pitch angle dispersion and i//(t) deter- 
mines the time variation of the beam. If the energy is much larger 
than the lower cut-off energy, z » 1, the distribution is power 
law with index -y — 1, thus, the flux spectrum (~ zf) is power 
law with index y. In the opposite case, z« 1, the distribution is 
power law with index 6 — 1 . The low energy index 5 is chosen to 
be 10 (see, e.g.. IZharkova & Gordovskvvl l2005a). while for the 
high energy index two values, 3 and 7, are considered. f D is the 
normalisation coefficient, which is chosen so that the energy flux 
of the injecting electron beam, 



F(s = s min ) = F 



/*/ 



dfiz nf{s = s min ,z,/J-), (16) 
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The distribution function is calculated in the following range 
of energies: z min < z < Zmax, where z min =0.1 (or 1.2 keV) and 
Zmax = 100 (or 1.2 MeV). The boundary conditions on energy 
are 



df(j,S,Z = Zmin,^) 

dz 

df(T,S,Z = Zmax,;") 

dz 



0, 



= 0. 



(18) 
(19) 



The boundary conditions on pitch angle are dMcClementsI 

1990) 



df(T,s,z,fi = 1) 
df(r,s,z,fi = -1) 



= 0, 
= 0. 



(20) 
(21) 



2.3. Integral characteristics of the electron distribution in the 
beam 

In the following sections we will numerically solve Eq. dTOb and 
calculate the following quantities for the electron beam: beam 
density (in cm -3 ), 

, 1 

n b (t, s) = F Q ^ ^| J dz f d» ViA (s) f (r, s, z, , (22) 

differential particle flux spectrum (in erg" 1 cm" 2 s" 1 ), 

l 

T» (r, s, z) = ^ fdfizA(s)f (r, s, z, fx) , (23) 
2£ 2 J 



mean particle flux spectrum (in erg 1 cm 2 s ') dBrown et alj 
I2003T) . 

J ds J dyar x (s)A(s)z,f (t,s,z,/j) 

(T n ) (r, z) = ^ s ^—l _ , (24) 

2 'fn- l (s)ds 

Smin 

angle distribution (in arbitrary units), 

(IN h (T.fl) 



dfi 



J ds J dztr l {s)A{s) ^f(T,s,z,n), (25) 



and energy deposition (or heating function) of the beam (in 

erg cm" 3 s" 1 ), 



Zmax 1 

/ ( r . s ) = (s)A(s) J dz J dfi |- fizf (t, s, z, /u) , 

(26) 
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where dz/ds is the electron's ene rgy losses with depth, which 
can be estimated as (Emslie 1980) 



10 



x 10 



dz _ I dz\ (dz\ _ 1 2 e 
ds \ds) c \ds) r yuz n 



(27) 



where two terms represent the collisional and Ohmic energy 
losses respectively. Coefficient A (s) = Bq/B{s) takes into ac- 
count the variation of the magnetic tube cross-section. 

3. Summary approximation method 

Let us combine the relative derivatives and rewrite the Fokker- 
Planck equation ( fTOb in the following form 
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df 
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1 -fi 2 d 2 f 
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+ 
d 2 f 



(28) 



Eq. (128b is solve d numerically b y using the summary ap- 
proximation method ( S amars kiil200 1 1) . This method allows us to 
study time dependent Fokker- Planck equation and it is different 
from the one used bv lZharkova & Gordovskvvl d2005bl) to solve 
the stationary problem. According to the summary approxima- 
tion method the four-dimensional problem is reduced to a chain 
of three two-dimensional problems. This is done by considering 
the three-dimensional differential operator at the right hand side 
of Eq.[10]as a sum of one-dimensional operators, each acting on 
the distribution function separately during one third of the time 
step. On each time substep the distribution function is calculated 
implicitly, hence, the numerical scheme is 



f 



/r+jAr _ f r = AT( p sLsf r +1 AT 



/ r+AT -r +iAT = Ar(^+, 



>2„ 



L2,)f 



(29) 
(30) 

: Ar (31) 



where L a are the finite difference operators that approximate 
the first order differential operators d/da. If the coefficient <p a 
is positive then the right difference scheme is used, i.e. L a f = 
/yar+A<* _ ytrj /Aar, otherwise the left scheme is used, i.e. L a f — 

(V - f a - Aa ) /Aa. The L^f = (f^ - 2f» + f^) I Aa 2 is 
the central difference that approximates the second order deriva- 
tive d 2 f/d/j. 2 . The computational grid has 200 nodes in the s 
dimension, 50 nodes in the z dimension and 30 nodes in the // di- 
mension. The nodes are distributed logarithmically in the s and 
z dimensions and linearly in the fi dimension. 

Eq. d29l ) together with boundary conditions forms a set of 
linear equations, which, after being solved, gives f T+ -A T from 
known f T . Distribution function / t+ " At is then used in Eq. ( f30T > 

to obtain f T+jAT - Finally, from Eq. OTT l we obtain / t+At which 
is, in turn, used in Eq. d29l on the next step. 

Electric field s is calculated on each time step according to 
Eq. (fT2] > where the distribution function is taken from the previ- 
ous step. Thus, the numerical scheme is not fully implicit. This 
means that, in order to avoid numerical instability, the time step, 
At, must be shorter than some critical value At c . In practice, 
the time step was determined by the trial-and-error method. For 
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Fig. 2: Electron density (a) and self-induced electric field (b) pro- 
files, where t is time passed after the injection is "turned on". 
Collisions and electric field are taken into account. The beam 
parameters, see Eq. ( fTBT ), are y = 3, F lop = 10 10 erg ctrT 2 s _1 
and Afi = 0.2. 



example, for the energy flux 10 10 erg crrr 2 s~' the time step is 
1 .7 x 10~ 4 s. It was found that when the energy flux is increased, 
the time step need to be decreased proportionally to keep the 
numerical scheme stable. 



4. Stationary injection 

Let us first recapture the major simulation results for the case of a 
stationary (steady) injection which helps to understand the phys- 
ical process of energy losses by beam electrons during their pre- 
cipitation. In our previous paper Zhark ova et al.l d2009l) we stud- 
ied hard X-ray emission and polarisation produced by a steady 
beam while in the current paper we focus more on the compari- 
son of electron precipitation results obtained for different models 
of magnetic convergence and on the energy deposition at various 
depths and times by electron beams with different parameters. 

The electron injection starts at t — and the simulation con- 
tinues until the stationary state is reached. The initial spectral 
index of the beam is chosen to be y — 3, the energy flux at the 
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Fig. 3: Beam electron density in different energy bands, if only 
collisions are taken into account. Beam parameters are the same 
as in Fig. [2] 



top boundary is F top = 10 10 erg cm 2 s 1 and the initial angle 
dispersion is A/i = 0.2. 

4.1. Relaxation time to a steady injection 

Let us first show how the system relaxes to the stationary state. 
Fig- H] shows the profiles of the electric field and beam density 
at different times. It is seen that the electric field relaxes some- 
what faster than the density. The relaxation time t r , after which 
the system becomes stationary, can be estimated as ~ 0.07 s. 
As it was stated in Sec. 12. li the self-induced electric field is cal- 
culated under the assumption that the beam current changes on 
time scales longer than the plasma collisional time. Since the 
collisional time in the corona is about 0.05 s and it decreases 
essentially with depth we may conclude that this assumption is 
correct. 

It was found that this relaxation time is longer for a stronger 
beam, which is the result of a smaller density of the ambi- 
ent plasma (Fig. []]) and, hence, a larger linear depth for the 
same column depth. For example, the depth 10 20 cm" 2 cor- 
responds to 2 X 10 8 cm and 11 x 10 8 cm for the beam en- 
ergy fluxes 10 10 erg cm" 2 s" 1 and 10 12 erg cm" 2 s" 1 respec- 
tively. This leads to a longer relaxation time for the atmosphere 
preheated by a stronger beam, e.g. for the beam energy flux of 
10 12 erg ctrT 2 s" 1 the relaxation time is found to be ~ 0.2 s. 

4.2. Effects of collisions and electric field (a B = 0) 

The local maximum, which appears on the density profile at 
the depth of about ~ 2 x 10 19 cm" 2 (Fig. |2ab . is caused by the 
beam deceleration while the flux of electrons remains nearly 
constant. After this depth most of the electrons leave the dis- 
tribution (thermalise) by reducing their energy below the en- 
ergy z m i n , and the density rapidly decreases at the depth of about 
~ 5 x 10 19 ctrT 2 (Fig.[2aJ». This corresponds to the stopping depth 
for the electrons with energies close to the lower cut-off energy 
(12 keV). The electrons with higher energies can travel deeper 
(see Fig. |3). In particular, it can be seen that electrons with en- 
ergies > 500 keV can travel down to the photosphere almost 



wi thout any energy losses. F ig. [3] confirms the results obtained 
bv lZharkova & Gordovskvvl d2006l) (see their Tab. 1). 

The depth profile of the electric field (Fig.l2bb has some no- 
ticeable variations at the column depths of about 4 x 10 18 cm" 2 , 
2 x 10 19 cm" 2 and 6 x 10 19 cm" 2 . These features are located at 
the same depths where the temperature profile (Fig. [TJ has sharp 
decreases leading to the sharp decreases of the plasma conduc- 
tivity given by Eq. (fT~3T >. However, the electron density fluctua- 
tion associated (through the Gauss's Law) with such an electric 
field profile is found to be of the order of 10" 4 cm" 3 and can be 
neglected. 

In Fig. [4] we plot the differential flux spectra at different 
depths for the forward (yu > 0) and backward (yu < 0) moving 
electrons. It can be seen that the self-induced electric field does 
not change the spectra of the downward (/j > 0) moving elec- 
trons but essentially affects the spectra of the upward (yU < 0) 
moving ones. Since the pitch angle diffusion due to the electric 
field is more effective for the lower energy electrons, the spec- 
tra of the returned electrons is enhanced at low and mid energy 
(Fig. a. 

The number of electrons returning back to the source plotted 
in Fig.|5a]is smaller compared to the case when the self-induced 
electric field is not taken into account. However, even without 
the electric field s, a number of electrons with fj. < is essential 
owing to the pitch angle scattering (second and third terms on 
the right hand side of Eq. (fT0b~). 

The heating function plot (Fig. l5bb shows that if a self- 
induced electric field is taken into account the heating is most 
effective at the upper column depth (2 x 10 19 cm" 2 ) compared to 
a much lower one (10 20 cm" 2 ) in the pure collisional b eam re- 
laxatio n, which is consistent with the results obtained by Emslie 
(1980). Indeed, the inclusio n of the electric fiel d decreases the 
stopping depth (IZharkova & Gordovskvvl 12006) and increases 
the number of returning electrons, thus, reducing the number of 
electrons at larger depths. All these factors lead to the upward 
shift of the heating funct ion maximum. The theoretical heating 
curve for pure collisions (Svrovatskii & Shmeleva 1972) is plot- 
ted in Fig. |5b] 

4.3. Effects of a magnetic convergence 

The converging magnetic field acts as a magnetic mirror and can 
essentially increase the number of the electrons that move up- 
wards. Let us determine how large the magnetic convergence 
parameter, should be to have any noticeable effect on the 
distribution of beam electrons. In order to do so we compare 
the terms in front of df/dfi in Eq. d28l . Since the collisional 
pitch angle diffusion is much smaller than the one caused by the 
electric field (see Fig. [5a]), we compare the effects of the mag- 
netic convergence with those caused by the electric field. The 
magnetic convergence effects are stronger if a B > 2s/(nz). To 
estimate this expression we plot the dimensionless ratio s/n in 
Fig. [6] The minimal value of the ratio s/n in the interval from 
s = Smin to the stopping depth of low energy (z = 1) electrons, 
~ 5 x 10 19 cm" 2 , is about 0.2. Thus, the magnetic convergence 
would be more effective than the electric field for the electrons 
with energies higher than the cut-off energy if a B > 0.4. High 
energy electrons can travel much deeper into the chromosphere 
(see Fig. [3]), where the ratio s/n can be as low as 2 x 10" 4 , thus 
the magnetic convergence would be more effective for them if 
a B > 4 x 10" 6 . 

In the following subsections we present the results of sim- 
ulations for different models of the converging magnetic field. 
These results are illustrated by the mean flux spectra plots 
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Fig. 4: Differential flux spectra of the beam electrons integrated over the positive and negative pitch angles. Beam parameters are 
the same as in Fig. [2] 



(Fig. [7aJ» for the upward (ju < 0) moving electrons, while the 
spectra of the downward moving electrons are found to be very 
close for all convergence models. The energy deposition profiles 
for different magnetic field approximations are shown in Fig.lTbl 

4.3.1. Exponential approximation of a magnetic field 
convergence 

Follow ing the approximation proposed by iLeach & Petrosiar] 
(1981), let us assume that the convergence parameter does not 
depend on depth 

<*b = fffio = const, (32) 



then the magnetic field variation is 



B(s) = B Q exp (a B0 (s - s min )) . 



(33) 



Suppose that the magnetic field at the depth s max is 1000 
times stronger than at the depth s m ; n , then ag ~ ln(1000)/i m ax = 
7.5 Xl0" 3 . As we discussed earlier the effect of the magnetic con- 
vergence with such low ag would be noticeable only for high 
energy (> 100 keV) electrons. This is clearly illustrated by a 
comparison of the mean flux spectra of the moving upward elec- 
trons and heating function obtained with and without magnetic 



convergence (see solid and dashed curves in Fig. |7). While the 
electric field returns mostly the low and mid energy electrons and 
makes the spectrum of the returning electrons softer (in compari- 
son with the purely collisional case), the magnetic mirror returns 
back the high energy electrons and makes their spectrum harder 
and similar to the initial power law. On the other hand, magnetic 
convergence reduces the heating at the deeper layers (Fig. l7bt . 
where it is caused by high energy electrons, because they were 
mirrored back to the corona. 



4.3.2. Parabolic approximation 

iMcClementsI dj992) suggested the following profile of a mag- 
netic field variation: 



B{s) = Bo 



1 + 



(^ l ^min) 



(34) 



If B(w)/5(s min ) = 1000 then s * W VflOOO) = 31.6. The 
convergence parameter is 



aiBis) 



2(s ^min) 
■^0 ^ ^min)^ 



(35) 
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Fig. 5: Pitch angle distribution (a) and heating function of the 
beam (b). Beam parameters are the same as in Fig. [2] and the 
magnetic convergence (for dotted curve) is given by Eq. i 



The magnetic convergence parameter ag at maximum is 1 /so ~ 
3.16 x 1CT 2 . This value is not high enough to affect all the beam 
electrons, but as is seen in Fig. [7] this convergence model, in a 
similar way as the previous one, increases the number of high 
energy electrons with negative fi an d reduces the heatin g of deep 
atmosphere layers. Note, also, that iMcClementsI (1 1 9921) consid- 
ered a constant plasma density, and in their model the profile of 
the magnetic field given by Eq. (l34l can be more effective, while 
in our case plasma density exponentially increases with depth. 
Here we assumed that the magnetic field changes according to 
Eq. < [34l > in the whole range of column depths. However it would 
be more appropriate to assume that the magnetic field variation 
is different in the corona and chromosphere. Such an approach 
is discussed in the following section. 



4.3.3. Hybrid approximation of magnetic field 

In this model we propose that is close to constant at small 
depth (in the corona) and tends to zero after some depth so (in 




column depth, cm 



Fig. 6: The ratio of the dimensionless electric field to the ambient 
plasma density as a function of column depth. Collisions and 
electric field are taken into account. Beam parameters are the 
same as in Fig. [2] 



the chromosphere): 

of B (s) = a B o 



Sq (^ ^min) 2 

then the magnetic field variation is 



(36) 



B(s) =Bq exp 



s 

J" a(s')ds' 



(37) 



Z?o exp a B oso arctan 



V so 



At shallow depths, where s <sc so, the magnetic field varies as 
B x Bq exp(a£o(s - s m ; n )), and at deeper layers, s » so, the 
magnetic field is constant, B « Bq exp(aBoSon/2). In most of the 
simulations (where it is not stated explicitly) we accept aso = 
10 and so = 0.2 (or 2.2 x 10 19 crrT 2 , which corresponds to the 
transitional region), this makes the ratio B(s max )/B(s m i n ) to be 
equal 23.1. 

Electrons with velocities inside the loss-cone are not re- 
flected by the magnetic mirror and reach the deep layers. For 
the current convergence model the critical pitch angle cosine of 
the loss-cone is n\ c = VI - B( s m i n ) / B( s max ) x 0.98. This means 
that for the accepted initial angle dispersion of 0.2 about 90% of 
the electrons are reflected back. 

As it is seen in Fig. [5b] the effect of magnetic convergence 
on the pitch angle distribution is similar to the effect of the 
electric field. However, since the electric field is most effec- 
tive for the electrons with fi = ± 1, the pitch angle distribution 
has a maximum at fi = 1 when the convergence is not taken 
into account. On the contrary, the magnetic field does not af- 
fect electrons moving along the field lines, thus, the angle dis- 
tribution of the upward moving electrons has a maximum at 
U m ~ -0.8 (Fig. |5a|). which is co nsistent with the conclusions 



Zharkova & Gordovskyy (2006). 



of 

As the convergence parameter a% is relatively high in this 
model, the whole spectrum of electron energy is affected (see 
Fig. [7aJ». The energy deposition profile for this magnetic field 
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Fig. 7: Mean flux spectra (a) of the upward (// < 0) propagating 
electrons and energy deposition (b) without (crosses) and with 
magnetic convergence given by Eq. ( f33b (solid curve), Eq. ( [34-b 
(dashed curve), Eq. < [37T > (dotted curve) and Eq. ( |38l (dot-dashed 
curve). Beam parameters are the same as in Fig. [2] 



approximation (Fig.lTbl indicates that the heating is only about 
30% of the heating produced in the case of constant magnetic 
field, which is because many of electrons are reflected by the 
magnetic mirror before they get into dense plasma. 

The profiles of electron density with different energies are 
plotted in Fig. [8] If only magnetic convergence is taken into 
account (Fig. [Sat , it can be seen that magnetic mirroring does 
not depend on the electron energy. Electrons the pitch angles 
of which are outside the loss-cone are turned back at depth 
~ 10 19 - 10 20 crrT 2 . The remaining electrons (the pitch angles 
of which are inside the loss-cone) can travel down to the lower 
boundary in the atmosphere. When the collisions are taken into 
account, electrons, especially those with low energies, lose their 
energy due to collisions (Fig. [8bl , It is important to compare 
collisional beam relaxation with and without magnetic conver- 
gence plotted in Figs.[8bland[3lrespectivelv. It is noticeable that 
the combination of the effects of collisions and convergence is 
stronger than the sum of the two separate effects. This occurs 



Fig. 8: Beam electron density in different energy bands. The 
magnetic convergence parameter is given by Eq. d36i . Beam pa- 
rameters are the same as in Fig. [2] 



because electrons with the initial pitch angles inside the loss- 
cone are scattered by collisions to pitch angles which fall out 
from the loss-cone, therefore, more electrons are returned back 
by the magnetic mirror. 

In Fig.[9]the results of simulations are presented for different 
magnitudes for the parameters aso and so of the magnetic con- 
vergence model given by Eq. ( 136*1 ). The increase of aso clearly 
affects the beam electrons by reducing the depth of their pene- 
tration and by increasing the number of returning electrons. It 
is obvious that the system would be sensitive to the variation of 
so if it is smaller than the penetration (stopping) depth. This is 
proven in Fig. [9a] 

4.3.4. Magnetic field model fitted to the observations 

Although the magnetic field cannot be directly measured in 
the solar atmosphere, there are some indirect techniques which 
allow the magnitude of the magnetic field to be estimated. 
The coronal magnetic field can be deter mined from radio 
observations of gyro-resonance emission (lLang et all 11993b 
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Fig. 9: Beam density (a) and energy deposition (b) as a function 
of column depth. The magnetic convergence parameter is given 
by Eq. ( T36T ) with various a bo and so- Beam parameters are the 
same as in Fig. [2] 



iBrosius & White! 120061) . In particular, IBrosius & Whitel d2006) 
suggest that the magnetic scale heigh t above sunspots, L g col = 
B/AB, is ~ 7 Mm. On the other hand. lKontar et al.l(l2008l) deter- 
mined the chromospheric magnetic field by measuring the sizes 
and heights of hard X-ray sources in different energy bands. 
They found that the chromospheric magnetic scale height, L^hr, 
is ~ 0.3 Mm. Assuming that the magnetic scale height, Lb, 
changes linearly with depth I, the convergence parameter is 



«bO) 



a bo 



(38) 



where aso = Ei/(ne nola AL^cor) = 15.7, and the dimension- 
less magnetic scale height as a function of the column depth is 



1-1 



(39) 



where the linear depth is l(s) oc J n(s)~ l ds. Since the magnetic 
field is defined as function of the linear depth the model depends 
on the density profile n(s) of the background plasma. 



Fig. 10: Electron density at impulsive injection as a function of 
column depth (a) and pitch angle (b). Collisions and the electric 
field are taken into account. The beam parameters, see Eq. < fT~5T >. 



are y — 3, F, 



top 



10 10 erg cm 2 s 



Afi = 0.2 and injection 



time, St= 1.7 x 10~ J s. 



As it is seen in Fig. [7a] such the magnetic convergence affects 
the electrons with all energies. This leads to about 20% reduction 
of the heating produced by beam electrons in comparison with 
constant magnetic field profile (Fig.lTbl. 



5. Impulsive injection 

Now we can consider the precipitation of electron beam inj ected 
as a short impulse. The early SMM dKiplinger et aill 19831) and 
recent CORONAS/IRIS dCharikov et al.ll2004l) observations re- 
veal millisecond impulses in the hard X-ray emission from solar 
flares. Furthermore, the electron accelerati on time in a recon- 
nectin g current sheet was recently shown ISiverskv & Zharkoval 
(2009) to be as short as 10 5 s. These facts suggest that the time 
scale of a beam of accelerated electrons may be rather short. 
In this section we study the evolution of such short impulse in 
the solar atmosphere. The injection time, 5t, is chosen to be 
1.7 x 10~ 3 s, which is much shorter than the relaxation time 
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Fig. 11: Mean flux spectra of the electrons injected as a short impulse. The beam parameters are the same as in Fig.fTOl 



t T ~ 0.07 s, that was found in Sec. 14. 11 (see Fig . [2j> . The default 
parameters of the beam are similar to the case of the stationary 
injection: the initial spectral index of the beam is y — 3, the max- 
imal energy flux at the top boundary is F top = 10 10 erg ctrT 2 s _I 
and the initial angle dispersion is A/j. = 0.2. Also the energy de- 
position profiles produced both by a softer beam (y = 7) and by 
a stronger beam (F top = 10 12 erg crrT 2 s _1 ) are obtained. 

The impulse injection, obviously, leads to a smaller density 
of electrons at a given depth in comparison with the stationary 
injection (see Figs. I2aland ll0al ). A smaller density results in a 
lower self-induced electric field. Thus, in the case of a short im- 
pulsive injection the electric field does not affect so much the dis- 
tributions. As a result, the only mechanism that can essentially 
increase the number of returning electrons is a magnetic con- 
vergence. As it was mentioned above in the current study we as- 
sume that the beam current is always compensated by the plasma 
return current, thus the self-induced e lectric field develop s im- 
mediately. However, as was shown by Ivan den Oordl ( 1 19901) . the 
neutralisation time of the beam current is of the order of the col- 
lisional time. Thus, the effect of electric field for short impulses 
can be even smaller than our estimations. 

Anisotropic scattering of beam electrons in collisions with 
the ambient plasma makes the pitch angle distribution more flat 
with time (see Fig. llObt . The electrons propagating downward 



reach depths with high density of the ambient plasma, lose their 
energy due to collisions and leave the distribution (become ther- 
malised) when their energy is less than z m j n . On the contrary, the 
returning electrons move into less dense plasma almost with- 
out losing any energy, but gaining it in the self-induced electric 
field. Thus, after some time the number of upward moving elec- 
trons can exceed the number of downward moving ones, which 
is clearly seen in Fig. llObl The angle distributions show that af- 
ter ~ 3.4 x 10 -2 s most of the downward propagating electrons 
are gone and the majority of electrons have [i < 0, i.e. they move 
back to the source in the corona. 

5.1. Energy spectra 

Since the first term at the right hand side of Eq. (I lot , which is re- 
sponsible for the energy losses due to collisions, is proportional 
to z~'^ 2 (where z is the dimensionless energy), one might expect 
that electron spectra would become harder with time. However, 
the downward moving electrons with higher energy reach the 
dense plasma faster and lose their energy faster than lower en- 
ergy electrons, which makes the energy spectra softer with time 
(Fig. |1 1 ab - The same is valid for the spectra of the upward mov- 
ing electrons. In this case, the high energy electrons escape the 
distribution faster by reaching the top boundary (s = Smm)- Due 
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Fig. 12: Energy deposition by a beam impulse. The beam parameters are the same as in Fig. [10] and the magnetic convergence is 
given by Eq. ( l3"6l l. 



to this effect the power law index can increase from the initial 
value 3 up to 4 during the beam evolution (Fig. II lab . 

Magnetic field on its own cannot change the energy of elec- 
trons. However, a converging magnetic field acts as a magnetic 
mirror and returns the essential part of electrons back to the 
source. As was shown above in Fig.[7a]for the high energy elec- 
trons the magnetic convergence is more effective than the elec- 
tric field and pitch angle diffusion. Thus, high energy electrons 
can quickly escape through the s = s m ; n boundary and the power 
law index can reach higher values than in the case with the con- 
stant magnetic field. For example, for the magnetic field profile 
given by Eq. d37b the power law index increases from the initial 
value 3 up to 8 (Fig. II let . If the convergence parameter is de- 
fined by Eq. d38l , the initial power law distribution converts to 
some kind of quasi-thermal distribution with an essential drop in 
high energies (Fig. 1 1 ldb . 

5.2. Energy deposition 

Fig. [12] shows the evolution of energy deposition, or heating 
functions, when different precipitation effects are taken into ac- 
count. In the purely collisional case (Fig. ll2al i the heating maxi- 
mum appears at the bottom boundary, moves upwards with time 



and vanishes near column depth ~ 10 iy - 10 20 cm . This evo- 
lution is consistent with stopping depths obtained for electrons 
with different energies (Fig. [3}. Indeed, the high energy electrons 
are the first to reach depths where the density is high enough to 
thermalise them. Less energetic electrons travel longer because 
of lower velocity but lose their energy higher in the atmosphere. 
Thus, the heating function maximum moves with the precipita- 
tion time from the photosphere (10 ms) to the lower (17 ms) and 
than upper (22 ms) chromosphere towards the stopping depth of 
the low energy electrons, where it vanishes. 

In a presence of the self-induced electric field (Fig. I12bb the 
heating by collisions becomes smaller than in the purely colli- 
sional case because some electrons are reflected by the electric 
field and do not reach dense plasma and, thus, do not deposit 
their energy into the ambient plasma. On the other hand, there is 
a second maximum on the heating function, which is also caused 
by the energy losses due to Ohmic heating by the self-induced 
electric field. This maximum does not move but grows in time 
with more electrons coming to the region with a high electric 
field (see Fig.l2bl. 

As it was shown in Sec. |4.3.3l for the magnetic convergence 
given by Eq. (|37] | only about 10% of the electrons can escape 
through the loss-cone and heat the deep layers. Thus magnetic 
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Fig. 13: Energy deposition by the impulse of a beam with y = 7. 
The other beam parameters are the same as in Fig.fTOl 



convergence substantially reduces the energy deposition at lower 
atmospheric levels due to the mirroring of the electrons back 
to the top, which shifts the heating maximum upwards to the 
corona (Figs. lT2"cl and ll2dY 

The heating function of a beam with the initial spectral index 
y — 7 is shown in Fig. [T3] Contrary to the y — 3 beam, the 
heating peak appears at shallower depths and moves downwards 
with time. Apparently, this occurs because of the number of high 
energy electrons being extremely low for a softer beam (y = 7), 
and, thus, the heating that they can produce at greater depths is 
too low to be noticeable. Also, the heating profile is narrower 
but higher and its maximum located higher in the atmosphere 
in comparison with y — 3 case. When the electric field is taken 
into account (Fig. |13bt the heating becomes stronger at smaller 
depths in the corona, in comparison to the pure collisional case 
(Fig. |13at . where it has maximum in the upper chromosphere. 



A more powerful beam, with energy flux 10 12 erg ctrT 2 s _1 , 
obviously deposits more of its energy in the ambient atmosphere 
(Fig. IT4T > than a beam with energy flux 10 I() erg crrT 2 s _1 . Two 
maxima are clearly seen on the heating function profile - one 
in the chromosphere, another one in the corona. If the mag- 



netic convergence is absent, the chromospheric heating is much 
stronger (Fig. I Mat . On the other hand, if the convergence is 
taken into account, then only about 10% of electrons can reach 
the chromosphere. Thus, the heating under the transition region 
is reduced by an order of magnitude, while the coronal heating 
remains nearly the same as in the case of the constant magnetic 
field (Fig. |14cb - Note, that in this case a different hydro-dynamic 
model is used to estimate the density and temperature of the 
ambient plasma. The discontinuity at the depth of 10 20 cm" 2 
(Fig. [141 is caused by a sharp increase of the ambient plasma 
density (see Fig.[TJ, which apparently corresponds to the transi- 
tion region. 

It can be noted that the evolution time of an electron impulse 
is longer for a stronger beam. This is a result of a smaller density 
of the ambient plasma, which leads to a longer relaxation time 
as was shown in Sec. 14. 11 



5.3. Bursts of hard X-ray emission 

In order to make a comparison with observations we calculate 
the intensity of hard X-ray emission produced by the injec- 
tion of a short electron beam. The bre msstrahlung cross-section s 
are taken in the relativistic form (see Gluckstem & Hullll 1953b . 
Fig. ll5al shows the time profile of hard X-ray intensity produced 
by the impulse with length 6t e = 1.7 ms. The timescale of the 
hard X-ray impulse, dt^, is about 20 ms, which is determined 
by the relaxation time of the atmosphere t r m 70 ms established 
earlier i n Sect. | 4.1| This times cale agrees closely with the obser- 
vations (Chari kov et alj |2004). Further simulations show that as 
long as St e <K tffhxr, the hard X-ray timescale depends only on 
the atmosphere parameters and does not depend on the length of 
the initial electron impulse. 

Evolution of the spatial profile of the hard X-ray intensity 
(Fig. ll5bt resembles that of the energy deposition (see Fig. I12bt . 
The emission starts at the bottom when high energy electrons 
reach this depth and gradually moves upwards. After reaching 



the depth ~ 2 x 10 cm the intensity the emission decreases 
and finally the emission vanishes. 



6. Conclusions 

By solving numerically the time-dependent Fokker-Planck equa- 
tion one is able to study the temporal evolution of the electron 
beam precipitation in the solar atmosphere and evaluation of 
the relaxation time required for the beam to reach the station- 
ary regime. For the beam with energy flux 10 10 erg cmr 2 s _1 
this relaxation time is 
factor of about 3 (or - 



- 0.07 s and it becomes longer by a 
0.2 s) for the beam with energy flux 

10 12 erg cirT 2 s _1 . 

The effect of the self-induced electric field during the station- 
ary beam injection is similar to that found in pre vious studies by 
Emslie ( 1980); Zha rkova & Gordovskvvl(l20 06). In particular, if 
the electric field is taken into account, then the maximum of the 
energy deposition profile is shifted upwards making the coronal 
heating stronger and the chromospheric heating weaker than in 
the case of pure collisional precipitation. 

We considered different models of a converging magnetic 
field to study the effectiveness of the beam electron refraction 
by a mag netic mirror. Magnetic field approximations used ear- 
lier bv lLeach & Petrosianl dl98lb : lMcClementsl dl992l) have the 
same spatial dependence in the corona and chromosphere. Even 
if the magnetic field in the photosphere is accepted to be 3 or- 
ders of magnitude higher than in the corona, such magnetic pro- 
files are shown to affect only the high energy electrons of the 
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beam. We propose a model where the magnetic field increases 
exponentially with depth in the corona and becomes constant in 
the lower chromosphere. Such magnetic field variation can affect 
the whole energy spectrum of electrons, while the ratio of photo- 
spheric/coronal magnetic field is as low as 23. Since the converg- 
ing magnetic field returns many electrons back to the source, the 
heating due to collisions and electric field is reduced by 70% in 
comparison with the constant magnetic field. We also considered 
the model based on indirect measurements of the magnetic field 
in the solar atmosphere. Such magnetic field variation can also 
affect electrons of all energies and reduce the collisional heating 
by 20% in comparison with the constant magnetic field profile. 

The further study is dedicated to the impulsive injections of 
electrons. In the simulation of impulsive injection the length of 
the impulse is chosen to be 1.7 x 10~ 3 s, which is much shorter 
than the relaxation time. It was found that the effect of the elec- 
tric field is considerably smaller for the short impulse than for 
the steady injection. 
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Initial energy spectrum of the injected impulse was power 
law. It was shown that during the evolution of the impulse the 
power law index increases in time. For example, if the joint ef- 
fects of the collisions and magnetic convergence are taken into 
account, the initial power law index of 3 can increase up to 8 and 
resembles some quasi-thermal distributions. 

The energy deposition profile is shown to depend on the ini- 
tial power law index. If the energy spectrum is hard (y = 3) the 
heating starts at the bottom end of the system due to the high 
energy electrons. In the case of soft (y = 7) impulse the num- 
ber of high energy particles is too low to produce any noticeable 
heating of the deep layers. On the other hand the higher layers 
are heated more effectively due to the higher number of the low 
energy electrons in the softer beam. 

We also compared the evolution of beams with different in- 
tensities. It was found that the difference in this case is mostly 
cased by the different density an d temperature profiles taken 
from the hydro-dynamic model (Zharkova & Zharkov 2007) 
(see Fig. [TJ. For example, the timescale of the impulse evolu- 
tion is longer the more intense the beam. 

If the timescale of the electron impulse is short enough, then 
the timescale of the hard X-ray emission is determined by the 
reaction of the atmosphere. Thus, it is of the same order of mag- 
nitude as the relaxation time, which, in turn, is of the order of 
10 ms and longer. This means that shorter electron impulses can 
not be detected by the hard X-ray observations. 
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